Quantum chiral phases in frustrated easy-plane spin chains 
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The phase diagram of antiferromagnetic spin-S chain with XY-type anisotropy and frustrating 
next-nearest-neighbor interaction is studied in the limit of large integer S with the help of a field- 
theoretical approach. It is shown that the existence of gapless and gapped chiral phases found in 
recent numerical studies [M.Kaburagi et al., J. Phys. Soc. Jpn. 68, 3185 (1999), T.Hikihara et 
al., J. Phys. Soc. Jpn. 69, 259 (2000)] is not specific for 5 = 1, but is rather a generic large-5* 
feature. Estimates for the corresponding transition boundaries are obtained, and a sketch of the 
typical phase diagram is presented. It is also shown that frustration stabilizes the Haldane phase 
against the variation of the anisotropy. 
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In recent few years, the problem of possible nontriv- 
ial ordering in frustrated quantum spin chains has at- 
tracted considerable attention 0-^]. Nersesyan et al. 
predicted [EJ that in anisotropic (easy-plane) antiferro- 
magnetic S = i chain with sufficiently strong frustrat- 
ing next-nearest-neighbor (NNN) coupling, a new phase 
with a broken parity appears, which is characterized by 
the nonzero value of chirality 



((S n x S n+ i) 2 



(1) 



note that the definition ([!]) differs from the other, so- 
called scalar chirality K oc S„_i • (S„ x S„+i) which 
is often discussed in the context of the isotropic spin 
chains ||. This prediction was made on the basis of 
the bosonization technique combined with a subsequent 
mean-field-type decoupling procedure. A similar conclu- 
sion can be reached Q] by means of a mean-field decou- 
pling of the quartic terms in the Jordan- Wigner trans- 
formed fermionic version of the Hamiltonian, in the spirit 
of the Haldane's treatment of spontaneously dimerized 
phase p|]. Up to now, however, this prediction for S = | 



was not confirmed in numeric studies [[l||2| . On the other 
hand, two different types of chiral ordered phases, gapped 
and gapless, were found numerically in 5 = 1 easy-plane 
frustrated chain |^|||. At present, to our knowledge, 
there is no theoretical analysis addressing the problem 
of chiral ordered phases for the 5 > 1 case. 

The aim of the present Letter is to study the generic 
large- S behavior of antiferromagnetic easy-plane integer- 
s' chain with frustrating NNN interaction, described by 
the following Hamiltonian: 

H = ^{(SnS n+ l)A+j(SA+2)A + D{S z n ) 2 } , (2) 
n 

Here (SiS 2 )a = SfS§ + 5f 5f + A5f5f , S„ denotes the 
spin-5 operator at the n-th site, the lattice spacing a has 
been set to unity, J > is the nearest-neighbor exchange 



constant, j > is the relative strength of the NNN cou- 
pling, and < A < 1 and D > are respectively the 
dipolar (inter-ion) and the single-ion anisotropies. 

We argue that the existence of gapless and gapped chi- 
ral phases found in [||,|| is not specific for 5 = 1, but 
is rather a generic large-5 feature. Estimates for the 
corresponding transition boundaries are obtained, and a 
sketch of the typical phase diagram is presented. As a 
side result, we also show that the domain of stability 
of the Haldane phase against the anisotropy variation 
grows when the frustrating coupling j is increased, in 
accordance with the numerical results P] . 

We use the well-known technique of spin coherent 
states [ ^2|Jl3| which effectively replaces spin operators by 
classical vectors (S+,5*) = 5(sin0 n e iv ", cos# n ) and in- 
corporates the quantum features by means of the path 
integral over all space-time configurations of (0,ip). The 
classical ground state of (||) is well-known: the spins al- 
ways lie in the easy plane (xy), i.e. 9 — S, for j < | the 
alignment of spins is antiferromagnetic, ip n — (p + nn, 
and for j > | a helical structure with incommensurate 
magnetic order develops, tp n — tpo ± (w — \o)n, where 
Ao = arccos(l/4j), and the ± signs above correspond to 
the two possible chiralities of the helix. 

In one dimension the long-range helical ordering is im- 
possible since it would imply a spontaneous breaking of 
the continuous in-plane symmetry; in contrast to that, 
the existence of the finite chirality K z n — (sin(<^ n +i — tp n )) 
is not prohibited by the Mermin- Wagner theorem. 

The classical isotropic system has for j > j three mass- 
less modes with wave vectors q = 0, q — ±<5, where 
6 = 7r — Ao is the pitch of the helix. The effective field 
theory for the isotropic case is the so-called £0(3) nonlin- 
ear sigma model, with the order parameter described by 
the local rotation matrix |l^Jll[]. The physics becomes 
simpler in presence of anisotropy since the modes with 
q = ±5 acquire a finite mass. Our starting point will be 
the following ansatz for the angular variables 9, ip: 



1 



<Pn 



n/2 +p n + (£ n e 
mi + ip n + (iu„e 



n5n 



+ w„e 



■iSny 2 
)/2- 



(3) 



We assume that the fluctuations p, £, u> are small and 
that they are smooth functions of ri, slowly varying over 
the characteristic distance Iq = 2ir/Xo; the same property 
is assumed for the function A„ = ipn+i — V J n which can 



be viewed as a dual variable to ij> 14 1. 

After passing to the continuum in the effective La- 
grangian L = J dxC = hSJ2 n (^ ~ cos6 n )d t tp n - H(9, ip) 
we average over Iq, making the oscillating terms disap- 
pear. The resulting expression for C is 



c = ns{p(d t m - iei 2 /4) + [w*(d t o + w(dt?)]/4} 

- JS 2 {V[X](l - |£| 2 /2) + A oP 2 + (A 1 /2)\w\ 2 ) 

- (JS 2 /4){M|£| 2 + FA|<U| 2 }, 

where the following notation has been used: 



(4) 



V[X] = jcos2X-cosX-U + (j/2)cos2X(d x X) 2 , 
Uq = j cos 2Ao — cos Ao, F = cos Ao — 4j cos 2Ao, 
A = D + A(l + j) - Uo, M = 2[D - (1 - A)U Q ], 
A\ = cos 2 A + j cos 2 2A - Uo . (5) 

Integrating out the "slave" fields p ~ —{%/2JSAo)d t ij), 
w ~ (h/2JSAi)dt£, and passing to the imaginary time 
y = ict, c = J S{2Aq) 1 / 2 /h, we obtain the effective Eu- 
clidean action 

A E = — 



90 
E 



l -Jj dxdy{\{d y ^f + V[\])(l- l -\tf 

d 2 X{(d,C)(d,O + m 2 \e}, (6) 



where (X 1 ,X 2 ) = (z.ic'f), d = JS{A 1 FA) 1 ' i /%, = 
d/dX^, and the constants go, E, mo are given by 



go 



E 



A FA 
At 



1/2 



M 
FA 



(7) 



Further, integrating out the massive £ yields the effective 
action for ip only, with a renormalized coupling g D g : 



9cs 



1 

9eS 

■■9o/ (I 



dx 



90^ 
2-kE 



dy{\{d y ^) 2 + V[\]) 
ln(l + A 2 /m 2 )), 



(8) 



where A = 7r is the lattice cutoff. A similar derivation 
may be carried out for j < 1/4: starting from the ansatz 
of the type m) with real £, w and S = it, one arrives at the 
same result (|), but with Ao set to zero in all quantities 
defined in (g), (@). 

We have mapped the initial quantum ID model to the 
2D classical XY helimagnet at effective "temperature" 
<7eff described by the effective action (g). The validity of 



this mapping is determined by the requirements go "C 1, 
\w\ <C |£| <C 1, p <C 1, which translate into 

5»(2A ) 1/2 , Ae"" i/9 ° < m < Sg /E (9) 

The first inequality above means that we are not allowed 
to consider large j > S 2 /2, and the second one requires 
the anisotropy to be within a certain range. We will 
be mainly interested in the behavior of for j close 
to the Lifshitz point 2l = \-> then the condition for the 
anisotropy transforms into 



C7r 2 ee" 27rS ^ < 3^/8 < 1 



(10) 



where /j, = 1 — A + 4D/3, e = \j — and the constant 
C = 1 for j < ji, and £ = 2 for j > , respectively. 

The model (g) possesses two basic types of topological 
defects Jjl| : (i) domain walls connecting regions of oppo- 
site chirality, and (ii) vortices, existing inside the domains 
with certain chirality and destroying the long-range heli- 
cal magnetic order (only quasi-long-range helical order is 
possible at finite g c s)- Thus one may expect two phase 
transitions: Ising-type transition ("freezing" of the do- 
main walls) which corresponds to the onset of the chiral 
order, and the Kosterlitz-Thouless (KT) transition (vor- 
tex unbinding) corresponding to the transition from the 
gapless chiral phase with algebraically decaying helical 
magnetic correlations (cos ip(x) cos ip(0)) oc (l/x) 9c!! ^ 2 ^ 
to the gapped chiral phase with only short-range he- 
lical order (but still with the long-range chiral order 
(k z (x)k z (0)) — > const, x — > oo). For the two transitions 
to be possible, one has to assume (later this assumption 
will be checked self-consistently) that the critical tem- 
perature of the Ising-type transition is higher that the 
corresponding temperature of the KT transition. 

Inside the phase with broken chiral symmetry one 
can set ip — ±Aoir + <j>±, then A ~ ±Ao + d x 4>± and 
V[X] ~ \F(d x (t>±) 2 . One obtains then the following esti- 
mate for the KT temperature: 



KT 



(n/2)VF 



(11) 



The equation g e s = g^ T determines the transition from 
the chiral gapless to the chiral gapped phase at j > jh , as 
well as the transition from the non-chiral (Ao = 0) gap- 
less XY phase to the non-chiral gapped Haldane phase 
a t j < 3l ■ Note that ( |ll|) is still valid away of the Lif- 
shitz point j = jl, since the field 4> remains smooth far 
from the vortex core, and the KT transition temperature 
is determined by the logarithmic divergence in the free 
vortex energy at large distances. 

In order to estimate the critical temperature of the 
Ising transition, let us first make some observations con- 
cerning the properties of chiral domain walls. The do- 
main wall (DW) energy can be easily calculated in the 
vicinity of the Lifshitz point, where A <C 1, so that the 
potential V[X] ~ (1/8){(A 2 - Ag) 2 + (^A) 2 } takes the 



2 



form of the ip 4 model, and one readily obtains the static 
DW solution A = Ao tanh {Ao(a; — xdw)} and the corre- 
sponding energy (per unit length in the y direction) 



E 



DW 



Ag/3, 



3 ~ h « 1 



(12) 



Further, it is easy to see that the chiral DW cannot 
move freely, since the infinitesimal displacement of the 
DW coordinate xdw would cause global change of the 
phase ijj at x — > oo. The DW can only "jump" by the in- 
teger multiples of tt/Xq, then the phase at infinity changes 
by the integer multiples of 2n. The jump by mr/Xo in- 
volves formation of n vortices bound on the DW, the 
elementary n = 1 jump is schematically shown in Fig. [I]. 
The energy per such a bound vortex can be estimated as 



Ebv 



Fln(7r/A ), A « 1 



(13) 
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FIG. 1. Elementary "kink" of the chiral domain wall in- 
terface, corresponding to the jump of the wall by tt/Xq. The 
arrows show the angle ip; the vortex in ip corresponds to the 
two "half- vortices" in the fields <f>± — ip =p Xqx living at the 
opposite sides of the boundary. The position of the domain 
wall is indicated with the dashed line. 



Ising 






disordered 










chiral 




/ \ 


gapped 


KT / 


\ 




XY \ / 


\ 


chiral 


gapless \ / 


KT 


gapless 



Jl J 

FIG. 2. Sketch of the phase diagram of the 2D helimagnet 
described by (^|). 

Since the Ising transition is governed by the discrete 
fluctuations of the DW interface, it is natural to use 



the so-called Muller-Hartmann-Zittartz, or the "solid-on- 
solid" approximation ]l6[] . In this approach the transi- 
tion temperature is determined by looking for the point 
where the free energy a of the DW interface becomes 
zero; a simple calculation yields the following equation 
for the critical coupling g = g].'. 

o~ = E DW - -7- In f 1 + d [cotanh^ - 1]} = , (14) 
do I 2g J 

where do — tt/{Xq^/F) is the characteristic size of the 
bound vortex in the y direction (in the derivation of 
( |l4| ) we have assumed that the distance along the y axis 
between two successive "jumps" should be greater than 
do). This equation can be solved numerically, and for 
j < 0.26 the solution is well fitted with the function 
g\ ~ 1.62Ao + 0.28A§, thus at j — ► the Ising transition 
temperature is larger than the KT one, gl > g^ T ~ §Ao, 
confirming the consistency of our assumption. 

Away from the Lifshitz point the above discussion of 
the Ising transition is no more valid, because the charac- 
teristic size of the bound vortex and the DW thickness be- 
come comparable with the lattice constant, and the con- 
tinuum description breaks down. However, it is known 
that E DW saturates at C DW w 0.87 for j > 0.8 @; 
one could also speculate that for j — > oo the energy of 
the bound vortex E^ v — > CbvVJ? where Ct v is some con- 
stant, and then from (Jl4|) one obtains g\. ~ CfwV^/lnj 
at j — > oo. On the other hand, according to ( JTl] ) 
9^ T ~ * ^VJ m the- same limit. Thus, one may expect 
that above a certain critical value of j the Ising transi- 
tion temperature g\ becomes lower than gf? T , and the 
gapped chiral phase disappears. 

The resulting conjectured phase diagram of the 2D he- 
limagnet (^J) is shown in Fig. |^ It should be mentioned 
that our picture of the transitions in 2D XY helimag- 
net disagrees strongly with that presented in ||l5| . In 
the latter work, using the arguments of fUjl , it was con- 
cluded that at low temperatures the vortices are bound 
by strings, which would suppress the KT transition and 
make the Ising transition to occur first with increasing 
the temperature. However, the argument of ||l7jl is ade- 
quate only for systems with broken in-plane symmetry, 
which is not the case here. Another point is that in the 
description used in [fl5[ the fields c/r^ , measuring the de- 
viations from the two different possible helix states with 
opposite chirality, are allowed to live and interact at the 
same space-time point, which, to our opinion, is rather 
unphysical. 

The above picture of the transitions in the 2D XY he- 
limagnet is now easily translated into the phase diagram 
of the frustrated spin chain, which is schematically shown 
in Fig. |^ for D = 0. Very close to jh, where too ^> A, 
which in terms of e = \ j— jt\ and fi = (1 — A) + |l? means 
e <C 3/Li/ (8£7r 2 ), the renormalization of the coupling con- 
stant is small, g c g ~ g$, and the transition boundaries 
are approximately given by 



3 



Kg 

ir 2 S 2 



D 



3 + 5A, 



(15) 



Here the coefficient Kc.c — 1 f° r the transition between 
gapless and gapped chiral phases, K C: h — 0.94 for the 
transition from the chiral gapped to the Haldane phase, 
and Kh-.xy — 2 for the Haldane- XY transition. One can 
see that the slope of the transition lines in the vicinity 
of ]l is very large (proportional to 5 2 ), and for large 5 
the boundaries mover closer and closer to the classical 
Lifshitz point j = ]l = j- 




chiral 
gapless 



FIG. 3. Schematic phase diagram of the anisotropic frus- 
trated spin chain (ph. 

At larger deviations from j — j^, when mo <C A, one 
has the following equations for the phase boundaries: 



M c = 



£n = 



2 5 2 



(16) 



which are valid for yfe — y/e^ 3> 1/5. One can see that 
the chiral gapped phase shrinks with increasing j. It is 
interesting to note that for j < jl the Haldane phase 
is stabilized by the frustration, in accordance with the 
numerical results |5). Further away from jl, when Ao 
becomes of the order of 1, the theory breaks down; how- 
ever, from the above arguments concerning the behavior 
of g\ we expect that the chiral gapped phase disappears 
above certain critical value of j. 

Certain limitations of the present theory should be 
mentioned. Our approach does not distinguish between 
integer and half-integer 5, since we have integrated out 
the out-of-plane components, and the only remaining 
topological charge, in-plane vorticity, plays no role. The 
topological term present in the full theory of the unit vec- 
tor field contains another quantum number, the so-called 
Pontryagin index; for j < jl this term is known [ fl3| to 
suppress the KT transition for half-integer 5, preventing 
the appearance of the Haldane phase. At j > jl there 
is no topological term |10 1 1 1 , and one may expect that 



state of a half-integer spin chain at sufficiently strong 
frustration is spontaneously dimerized ||, and our ap- 
proach does not allow one to capture this feature. An- 
other limitation is that we cannot describe the hidden 
(string) order in any way, and thus it is not possible to 
analyze the coexistence of the string order and chirality 
in the gapped chiral phase observed in Refs. |4||| nor 
to study the transition to the so-called double Haldane 
phase characterized by the absence of the string order 
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the KT transition for j > survives also for half-integer 
5. However, this point is not so clear since the ground 
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